Stark deceleration of CaF molecules in strong- and weak-field seeking states 
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We report the Stark deceleration of CaF molecules in the strong-field seeking ground state and 
in a weak-field seeking component of a rotationally-excited state. We use two types of decelerator, 
a conventional Stark decelerator for the weak-field seekers, and an alternating gradient decelerator 
for the strong-field seekers, and we compare their relative merits. We also consider the application 
of laser cooling to increase the phase-space density of decelerated molecules. 
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I. INTRODUCTION 

A Stark decelerator is a device for slowing down pulses 
of cold molecules [T] . The force that acts on the molecules 
is the spatial gradient of their Stark shift arising from the 
inhomogeneous electric field of the decelerator. The abil- 
ity of a Stark decelerator to deliver cold polar molecules 
with a centre-of-mass velocity at or near zero has proven 
exceedingly useful for a wide variety of experiments in- 
cluding high precision spectroscopy [2] , measurements of 
the lifetimes of long-lived molecular states "3", '4', low tem- 
perature collision studies [SJ IH] , and tests of fundamental 
physics [ZIIH]. 

To date, only relatively light polar molecules in weak- 
field seeking states have been decelerated to low speed by 
the Stark deceleration method [U [9Hl4] . For many ap- 
plications it would be useful to extend the range of slow 
molecules available by applying the method to heavier 
molecules. This is particularly relevant for molecule- 
based tests of fundamental physics, such as the mea- 
surement of the electron's electric dipole moment [15] 
and tests of parity violation in nuclei |16j and chiral 
molecules [T^, all of which require heavy molecules to 
reach high precision and could benefit from the use of 
slower molecules |18) . 

The kinetic energy to be removed by the decelerator is 
proportional to the molecular mass, so heavier molecules 
require more deceleration stages. Fortunately, the decel- 
erator focusses the bunches of molecules both longitudi- 
nally [il9 and transversally [2^, so the bunches do not 
spread out as more stages are added. This important fo- 
cussing property applies only to molecules in weak-field 
seeking states, since they are naturally attracted to the 
axis of the machine where the field is lowest. Focussing of 
strong-field seekers in all three directions is not possible 
because of the impossibility of creating an electric field 
maximum in free space |21) . Unfortunately, the low- lying 
states of heavy molecules are all strong-field seeking in 
the large electric fields needed for a decelerator. This is a 
major obstacle to extending Stark deceleration to heav- 
ier molecules. There are two possible ways to circumvent 
this difficulty. One option is to prepare the molecules in a 
higher-lying rotational state which does have a weak-field 
seeking component at the relevant fields, and then decel- 
erate them in the usual way. The heavier the molecule 



the greater the rotational excitation will need to be in or- 
der to find a suitable weak-field seeking state. The second 
option is to an use an alternating gradient (AG) decel- 
erator [22H24] which works with molecules in strong-field 
seeking states. In this type of decelerator each decelera- 
tion stage focusses the molecules in one transverse direc- 
tion and defocusses them in the other, the focussing and 
defocussing directions alternating from one stage to the 
next. This series of focussing and defocussing lenses can, 
when properly arranged, provide net focussing in both 
transverse directions. An AG decelerator can be used for 
molecules in any state, including the ground state which 
always has the largest Stark shift, and its application to 
heavy molecules has been demonstrated [23, 25 . 

In this paper we describe the deceleration of CaF 
molecules in both weak-field seeking and strong-field 
seeking states. For the former case, we use CaF molecules 
in a weak-field seeking component of the fourth rotation- 
ally excited state {N = 4) so as to maximize the Stark 
shift, and we use a decelerator that follows a conventional 
design - we call it a WF-decelerator. For decelerating 
strong-field seeking CaF we use molecules in the ground 
rotational state and an AG decelerator. We compare 
these two strategies for decelerating heavy molecules. We 
also consider the prospects of laser cooling CaF so as to 
increase the number and phase-space density of deceler- 
ated molecules. 



II. EXPERIMENTAL SET-UP 

A pulsed supersonic beam of CaF is produced by 
laser ablation of Ca into a pulsed supersonic carrier 
gas (Ar, Kr or Xe) containing a small fraction of SFg 
|26j . The translational and rotational temperatures of 
the molecules are both approximately 3 K. Their speed 
is usually a little higher than the terminal supersonic 
speed of the carrier gas. The molecules pass through a 
skimmer, then through the decelerator, and are finally 
detected 805 mm from the source by time-resolved cw 
laser-induced fluorescence on a chosen rotational compo- 
nent of the - X2S+(0 - 0) transition at 606 nm. 
The experiment runs at a repetition rate of 10 Hz. Fur- 
ther details of the methods we use to produce and detect 
cold CaF molecules are given in p7j . 
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FIG. 1; Photographs of the two decelerators. (a) Three stages 
of the AG decelerator. (b) The first few stages of the WF 
decelerator. 

The AG decelerator, shown in Fig.[l|a), consists of 21 
deceleration stages, each of which is also a lens that fo- 
cusses in one transverse direction and defocusses in the 
other. Each of these lenses is formed by a pair of stain- 
less steel cylindrical electrodes, 14 mm long and 6 mm in 
diameter, with hemispherically-rounded ends of radius 
3 mm, their axes parallel to the z-axis, and their surfaces 
separated by 2 mm. The centre-to-centre spacing of the 
lenses (along z) is 24 mm, and successive lenses are ro- 
tated through 90° about the z-axis. When the electrode 
pair lies in the a;z-plane the lens focusses strong-field 
seeking molecules along y and defocusses them along x. 
The following lens has its electrodes in the yz-plane and 
does the opposite. Each electrode is attached by a pair 
of dowels to one of four 16 mm-diameter stainless steel 
support rods, 504 mm in length. Each support rod is 
attached by insulating macor stand-offs to two stainless 
steel support rings, one at either end of the decelerator. 

The WF decelerator, shown in Fig.[l|^b), consists of 
100 deceleration stages. Each stage is formed by two 
parallel stainless steel cylindrical electrodes of diameter 
3 mm, with their axes parallel to either x or y, and their 
surfaces separated by 2 mm. The centre-to-centre spac- 
ing of the stages is 6 mm along z, and successive stages 
are rotated through 90° about the z-axis. When the elec- 



trodes have their axes along x they focus weak-field seek- 
ing molecules along y, and do nothing along x. For the 
electrodes aligned along y the opposite is true. Each elec- 
trode is pushed into one of four 16 mm-diameter stainless 
steel support rods, 594 mm in length, and these support 
rods are attached by insulating alumina stand-offs to two 
stainless steel support rings. 

For both decelerators, four independent 20 kV switches 
are used to switch the high voltages applied to each of the 
support rods between two values, ±Vhi and ±VlOj with 
rise and fall times of approximately 500 ns. The support 
rings are grounded. For the AG decelerator, ±Vlo is 
zero. To avoid nonadiabatic transitions in the WF decel- 
erator i28], Vlo is not zero but is still far smaller than 
Vhi- After high voltage conditioning, both decelerators 
are able to support ±20 kV across the 2 mm gaps between 
electrode surfaces. 



III. PROPERTIES OF THE TWO 
DECELERATORS 

As is usual, we introduce a reduced position 6 — -kz/ L 
(modulo 27r) where L is the distance between deceleration 
stages. The decelerators can be in one of four high volt- 
age configurations: (i) even stages at ±Vhi, odd stages at 
±Vlo, (ii) odd stages at ±Vhi, even stages at ±Vlo, (iii) 
all stages at ±Vlo, and (iv) all stages at ±Vhi- The WF 
decelerator is switched back and forth between configu- 
rations (i) and (ii) in a sequence generated such that a 
molecule with a chosen initial speed Vi reaches the same 
point relative to the periodic array every time the decel- 
erator is switched. The value of at this point is termed 
the synchronous phase angle, 4>q [T5]. The switching se- 
quence for the AG decelerator is slightly more compli- 
cated - the decelerator is switched into state (i) when the 
synchronous molecule reaches the position Zon inside an 
even stage, to state (iii) when it reaches the position Zoff , 
to state (ii) when at Zon in an odd stage, and back to 
state (iii) when it reaches Zoff again. 

In Fig. [2] we compare the basic properties of the two 
decelerators. Here, we take the electrodes to be charged 
to ±20 kV so that the field on the beamline is approx- 
imately 200kV/cm, and we take the CaF molecules to 
be in the {N, M) = (0, 0) state for the AG decelerator 
and in the (4, 0) state for the WF decelerator. Figure 
[2ja) shows the Stark potential as a function of position 
(z) along the beamline for configurations (i) and (ii) . Al- 
though the electric field on the beamline is approximately 
the same in the two decelerators, the depth of the Stark 
potential is about 4 times higher in the AG than in the 
WF because of the larger Stark shift of the ground state. 
However, the stages are packed 4 times closer together in 
the WF and so the deceleration per unit length of decel- 
erator is approximately the same in the two cases. Figure 
[2jb) shows the longitudinal acceptance areas of the two 
decelerators, calculated analytically in both cases using 
the relevant effective potentials [121 [30] ■ In both cases, 
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FIG. 2: Key information for deceleration of CaF in the two 
decelerators: WF (solid red lines) and AG (dashed blue lines). 
Both decelerators are charged to ±20 kV. The molecules are 
in the (A'', M) = (4, 0) rotational state for the WF decelera- 
tor and in the (0, 0) state for the AG decelerator. (a) Stark 
shift versus longitudinal position with even stages charged 
and odd stages grounded (bold) and vice versa (thin), (b) 
Longitudinal acceptance when the energy loss per stage is 
half the maximum possible value, (c) Transverse acceptance 
when the energy loss per stage is half the maximum possible 
value. For the AG decelerator we show the acceptance area 
at the entrance of a focussing lens, we took an effective lens 
length of 5.5 mm, an effective drift length of 18.5 mm and a 
forward speed of 430 m/s. 



the energy loss per stage is chosen to be half the maxi- 
mum possible value. The longitudinal acceptance of the 
AG is larger than the WF because the stages are longer 
and the the Stark shift is larger. Figure ^c) shows the 
acceptance areas of the two decelerators in one of the two 
transverse directions. For the WF decelerator this is ob- 
tained by calculating the mean transverse spring constant 



PU] , while for the AG decelerator it is found by using the 
Courant-Snyder formalism for AG focussing j24 l l29 t [30 ] . 
We find that the transverse acceptance areas are approx- 
imately the same for the two decelerators. In the AG de- 
celerator each lens focusses/defocusses strongly, because 
of the large Stark shift. The net effect of the alternating 
lenses is to focus the molecules, and this can be described 
as an effective focussing lens whose power is consider- 
ably smaller than that of the individual lenses In 
the WF decelerator the molecules are focussed through- 
out, but the focussing is relatively weak because of the 
smaller Stark shift and because, whenever the molecules 
approach a region of strong field where the focussing is 
strongest the decelerator is switched to make the field 
small again. 

We have used approximate methods to calculate the 
phase-space acceptance areas shown in Fig. [2] In partic- 
ular we have used an effective potential to describe the 
longitudinal motion of molecules relative to that of the 
synchronous molecule, we have assumed that the trans- 
verse forces are harmonic, and that the longitudinal and 
transverse motions are uncoupled. This last approxima- 
tion is particularly poor, the motions in the longitudinal 
and transverse directions being quite strongly coupled in 
both the WF decelerator [5^ and the AG decelerator [H] . 
The true acceptances will be smaller as a consequence, 
but still of the same order of magnitude as obtained from 
these approximate methods (see |18| for example.) 



IV. AG DECELERATION RESULTS 

In this section we discuss the deceleration of ground 
state CaF from an initial speed of Vi = 433 m/s. The 
AG decelerator is used with Vhi = 20 kV, Vlo = 0, 
= mm and Zog = 6 mm, as indicated in Fig.[2ja). 
The deceleration switching sequence is applied on every 
even-numbered pulse of the experiment, while on every 
odd-numbered pulse the decelerator is turned off. Figure 
shows how the time-of- flight profile of the molecules 
changes as the number of deceleration stages used in the 
experiment, n, is increased. It is the last n stages that 
are used, the decelerator being off until the synchronous 
molecule is n stages from the decelerator exit. To elimi- 
nate the effects of drifts in the source flux, each profile has 
been normalized to the amplitude of the corresponding 
decelerator-off profile obtained (almost) simultaneously. 
The zero of time is the arrival time of the synchronous 
molecule (with speed Vi) when the decelerator is off. Af- 
ter a few stages the phase-stable molecules have been 
bunched about the synchronous molecule, resulting in a 
narrow peak in the measured time-of- flight profiles. As n 
is increased the narrow peak moves to later arrival times, 
showing that the speed of this bunch has been reduced. 

Figure [sj^b) shows the time-of-flight profiles obtained 
by simulating the motion of molecules through the decel- 
erator, using a field map obtained from a finite element 
model and removing any molecule that crashes into an 
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FIG. 3: (a) Points: Measured time-of-flight profiles of ground- 
state CaF witli tlie last n stages of the AG decelerator being 
used. The parameters are Vm = 20 kV, Vlo ~ 0, 2:011 = mm, 
ZoS = 6 mm, Vi = 433 m/s. The bottom profile was taken with 
the decelerator turned off. Lines: Triple Gaussian fits to the 
data, (b) Points: Simulated time-of-flight profiles to match 
the experiments. Lines: Triple Gaussian fits. 



electrode. The simulation results agree fairly well with 
the experimental data. To make a quantitative com- 
parison, we fit each of the experimental and simulated 
datasets to a model which is a sum of three Gaussians. 
In this model, one Gaussian (Gi) represents the broad 
background of undecelerated molecules, a second (G2) 
represents the decelerated bunch, and a third (6*3), with 
negative amplitude, represents the hole in the undeceler- 
ated distribution where the decelerated bunch would oth- 
erwise have been. These fits are shown as lines in figure 
[3j We see that this model fits well to the data, except for 
some modulations in the undecelerated background that 
the model does not capture. These modulations, which 
are due to bunching of molecules in the neighbouring de- 
celeration stages, are most clearly seen in the simulated 
data for small n, but are also present to a lesser extent 
in the experimental data, and in the data for higher n. 
As one would expect, the central arrival times of Gi and 
G3 have no significant dependence on n, while the arrival 
time of G2 increases with n. This is shown in Fig.|4]^a) 
where we plot, as a function of n, the delay in arrival 
time between the decelerated bunch and the arrival time 
of the synchronous molecule when the decelerator is off. 
Note that for the relatively small reductions in speed ob- 
tained here, the narrow bunch of molecules represented 
by G2 is actually a mixture of phase-stable molecules 
and some molecules that are not phase-stable but have 
not yet separated from the phase-stable bunch. The ex- 
perimental delays are similar to the predictions of the 
simulations but deviate for n — 19 and 21 where the sim- 
ulations predict a larger delay than is measured. We do 
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FIG. 4: Comparison of experimental (filled circles, red) and 
simulated (open circles, blue) results for the decelerated 
bunch, as determined by the parameters of a three-Gaussian 
fit. The error bars are the uncertainties on the fitted param- 
eters; when not visible, they are smaller than the size of the 
points, (a) Delay in arrival time as a function of n. (b) Num- 
ber of molecules in the decelerated bunch as a function of n, 
normalized to the total number of molecules measured when 
the decelerator is off. 



not know the reason for this discrepancy. In the exper- 
iments for n = 21 the phase stable bunch is decelerated 
from 433 m/s to 399 m/s, corresponding to a removal of 
2.1 THz of kinetic energy, or 15% of the initial kinetic 
energy. 

Figure |4]Jb) shows how the area of the decelerated 
bunch (i.e. of G2) depends on n for both the experimen- 
tal and simulated data. We have normalized these areas 
to the total area of the signal obtained when the decel- 
erator is off. The simulations predict an increase in the 
signal ratio up to n = 13 and then a decrease at higher 
n. When n is small, the distance from the source to 
the effective entrance of the decelerator is large, and the 
range of transverse speeds that can enter the decelerator 
is smaller than the range the decelerator can accept. The 
decelerator's acceptance is thus not filled. As n increases 
the effective entrance moves closer to the source, the ac- 
ceptance is more completely filled and so the number of 
molecules in the decelerated bunch increases. At around 
n = 13 the acceptance is completely filled and further 
increasing n now reduces the signal ratio because there 
is a loss mechanism at work in the decelerator, as we dis- 
cuss below. The measured signal ratio for the decelerated 
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bunch is always smaller than the simulations predict, re- 
maining constant up to n = 15, and then dropping. This 
is the behaviour we would expect if the continuous losses 
within the decelerator are larger than expected. We note 
that there is a similar reduction with increasing n in the 
transmission of the undecelerated molecules. 

The transverse forces at the ends of each lens are re- 
sponsible for these continuous losses. In an ideal AG 
decelerator, where the transverse forces are linear in the 
off-axis displacement with equal and opposite force con- 
stants, the transverse acceptance does not change as more 
lenses are added, at least not once the decelerator length 
is longer than the wavelength of the macromotion. As 
has been discussed previously [24j, the ends of the AG 
decelerator electrodes are detrimental to proper AG fo- 
cussing because the magnitude of the defocussing force 
constant ky = —dFy/dy greatly exceeds the magnitude 
of the focussing force constant — —dF^jdx near these 
ends, Fx,y being the components of the force in the fo- 
cussing (x) and defocussing (?/) directions. This is shown 
in Fig.[5ja), where we plot the quantity —ikx + ^y) as a 
function of position along the axis, z. For reference, the 
electric field along z is shown on the same plot. While 
the focussing and defocussing force constants are well- 
balanced inside the lens, the defocussing is far stronger 
than the focussing in the region near the ends of the 
lens. Figure [sjb) demonstrates the impact of these end 
effects. Here, we have simulated AG deceleration with 
n = 17, Vhi = 20 kV, Vlo = 0, Zon = mm, ZoS = 6 mm 
and Vi = 436 m/s. The speed distribution entering the 
decelerator is a Gaussian centred at 436 m/s. From the 
simulation we can select out those molecules that are 
transmitted by the decelerator and then look at their 
distribution of initial speeds. It is this distribution that 
is plotted in Fig.[5];b). We see that the transmission de- 
pends very strongly on the speed. We draw attention to 
molecules with three specific velocities, 420 m/s, 430 m/s 
and 440 m/s, indicated by the vertical lines, where the 
transmission is bad, good and bad, respectively. Each 
of these molecules passes through a 6 mm- wide region of 
the first lens during the time when it is switched on. For 
each molecule, the centre of this region is indicated by a 
vertical line in Fig.jsj^a). The molecules that are poorly 
transmitted are the ones that see the two bad regions of 
this first lens. The simulations show that these molecules 
repeatedly experience the bad regions as they move from 
one lens to the next, and that is why they tend to be lost. 
The 430 m/s molecule sees only the good part of the first 
lens, and to a large extent avoids the bad regions of the 
other lenses too, and so is transmitted well by the decel- 
erator. Unfortunately, the bad region at the exit of each 
lens is also the region where the molecules are deceler- 
ated, and so the decelerated molecules necessarily have 
to move through this bad region and tend to be thrown 
out of the decelerator for this reason. 

Our simulations use an accurate map of the field, cal- 
culated numerically using a finite element model, but the 
experimental data shows more molecule loss than in the 
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FIG. 5: (a) The quantity —{kx + ky) plotted as a function of 
distance z along the AG decelerator, indicating the regions 
of the lens where the defocussing is far stronger than the fo- 
cussing. For reference, the electric field as a function of z is 
also plotted (on a different scale), (b) Simulation result show- 
ing the distribution of initial velocities that are successfully 
transmitted through the decelerator. Molecules with three 
velocities are indicated by vertical lines and are discussed in 
the text. Their central positions inside the first lens, during 
the time that it is turned on, are shown by the equivalent 
vertical lines in (a). 



simulations. We do not know the reason for this. Other 
AG deceleration and focussing experiments also report 
greater losses than expected, the suspected cause being 
electrode misalignments j22j [25, i31j ,32] . 

V. WF DECELERATION RESULTS 

In this section we discuss the WF deceleration of CaF 
molecules that emerge from the source in the {N, M) = 
(4,0) state and with an initial speed of 340 m/s. For 
these experiments, it is important to suppress nonadi- 
abatic transitions to other (4, M) states which may be 
driven by the rotation of the electric field when the de- 
celerator switches, especially in regions where the electric 
field is small |28| . These unwanted transitions can be 
suppressed by ensuring that the Stark splitting between 
the various M states is much larger than the angular 
frequency at which the field rotates. For this reason, in- 
stead of switching between high voltage and ground, we 
switch between Vm and a bias voltage of Vlo — 1 kV. 
This size of bias voltage eliminates losses due to nonadi- 
abatic transitions and makes no significant difference to 
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FIG. 6: (a) Measured time-of-flight profiles of CaF in tlie 
A*' = 4 rotational state following deceleration in the WF de- 
celerator. The zero of time is the arrival time of the syn- 
chronous molecule when the decelerator is off. For all the 
data, Vi = 340 m/s, Vlo = 1.0 kV. The high voltage (Vhi) 
and synchronous phase angles {(j)o) for the individual traces 
are (i) 18 kV, 15°, (ii) 18 kV, 30°, (iii) 18 kV, 45°, (iv) 18 kV, 
60°, (v) 20 kV, 60°, (vi) 20 kV, 75°. The inset shows the de- 
celerated bunch for trace (vi) in more detail, (b) Numerical 
simulations of these same experiments but with Vlo = and 
non-adiabatic transitions neglected. 



the energy loss per stage. 

Figure [6] shows the results of these experiments, along 
with the associated numerical simulations. The laser- 
induced fluorescence signal from molecules in iV = 4 is 
plotted versus arrival time for several different values of 
the synchronous phase angle and the high voltage Vhi. 
The zero of time is the arrival time of the synchronous 
molecule when the decelerator is turned off, or when it 
is operated at 0o = 0. As (f>Q increases the deceler- 
ated bunch moves to later arrival times, corresponding 
to slower speeds, but also reduces in amplitude because 
the longitudinal phase space acceptance of the deceler- 
ator decreases as 0o increases. In the uppermost trace, 
we have set (po = 75° and Vhi — 20 kV. In this case 
the molecules are decelerated to a final speed of 224 m/s, 
corresponding to an energy reduction of 4.9 THz or a re- 
moval of 57% of the initial kinetic energy. The simulated 
time-of flight profiles shown in Fig.[6][b) agree well with 
the experimental data. In particular, the measured ar- 
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FIG. 7: Results from the WF decelerator operated at Vki = 
18 kV. (a) Points: Measured kinetic energy loss for the de- 
celerated bunch versus synchronous phase angle. Line: Ex- 
pected energy loss calculated from the Stark potential, (b) 
Points: Measured number of decelerated molecules versus 
synchronous phase angle, normalized to the number at (j)o — 0. 
Line: Longitudinal acceptance versus synchronous phase an- 
gle. 



rival time and size of the decelerated bunch are the same 
as predicted by the simulations. 

Figure [t] (a) shows how the measured energy loss 
changes with when Vm = 18 kV. Here, we have deter- 
mined the energy loss from the measured arrival time of 
the decelerated bunch and the measured distances from 
source to decelerator and decelerator to detector, by as- 
suming that the decelerator applies a constant force to 
the molecules over its entire length. Since there are many 
deceleration stages, each applying the same net deceler- 
ation, this is a good approximation. The energy loss 
follows the prediction obtained from the Stark potential 
plotted in Fig.|2];a). In Fig.[7];b) we show the number of 
decelerated molecules as a function of 4>o, normalized to 
the number at (f)Q — and compare this to the (similarly 
normalized) longitudinal acceptance. As 00 increases the 
number of decelerated molecules decreases because the 
longitudinal acceptance decreases. While the experimen- 
tal results follow the trend of the longitudinal acceptance, 
there is some deviation from this simple model. This is 
because the transverse acceptance also has some depen- 
dence on (f>Q, through the coupling of the transverse and 
longitudinal motions. This coupling tends to reduce the 
number of decelerated molecules at low phase angles and 
to increase it at high phase angles [50], as we indeed ob- 
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serve in Fig.[7]jb). 



VI. PROSPECTS FOR COMBINING STARK 
DECELERATION AND LASER COOLING 

Laser cooling of a diatomic molecule has recently been 
demonstrated for the first time [33] ■ Laser cooling is fea- 
sible if the molecule has a short-lived state that can be 
excited with a convenient laser and that decays (exclu- 
sively) to the ground state with a Franck-Condon factor 
close to unity. Then, only a few excitation lasers are 
needed to address all relevant vibrational components of 
the electronic transition [M]- The A^Hi/a - X^I]+(0 - 0) 
transition of CaF has all the required properties for laser 
cooling [2T. It may be possible to both cool and decel- 
erate the CaF beam using a frequency-chirped counter- 
propagating laser. To decelerate from 340 m/s requires 
the molecule to scatter about 30000 photons, and it is not 
yet clear how many of the vibrational transitions would 
need to be addressed to achieve this since the Franck- 
Condon factors are not known well enough. An alterna- 
tive approach is to combine laser cooling and Stark de- 
celeration. Laser cooling prior to deceleration increases 
the phase space density of the molecular pulse and so in- 
creases the number of molecules that can be captured by 
the decelerator. The decelerator is then used to bring the 
molecules to rest. Since the cooling only needs to change 
the molecule speeds by about 10 m/s to compress the ve- 
locity distribution substantially, the required number of 
scattered photons is greatly reduced. 

To estimate how effective laser cooling would be in 
increasing the number of decelerated molecules, we sim- 
ulate a possible experiment. In this simulation, the laser 
light consists of two main components, one at 606 nm 
resonant with the P(l) line of A{v' = 0) - X{v" = 0), 
and the other at 628 nm resonant with the P(l) line of 
A{v' = 0) — X(w" — 1). Sidebands are applied to ad- 
dress the 4 hyperfine components of each of these tran- 
sitions, so that there are 8 discrete frequencies in total. 
A small magnetic field is used to prevent optical pump- 
ing into a dark state. The only remaining dark states 
are the vibrational states with v" > 1, and the branch- 
ing ratio for the v' = state to decay into these dark 
states is assumed to be approximately 10~^. Neglecting 
coherence, the internal dynamics of the molecule inter- 
acting with this set of laser beams is calculated using a 
rate model that includes all the Zeeman sub-levels (24 
in the ground state, 4 in the upper state). For the de- 
cay rate and the — Franck-Condon factor we use the 
measured values [27], r/(27r) = 8.3 MHz, and Z = 0.987. 
The branching ratios between all the sub-levels are calcu- 
lated using the angular momentum algebra written down 
in [57]. To include the influence of the magnetic field, 
which does not fit comfortably in a rate model, we add 
terms to the rate equations that damp out differences 
in population between Zeeman sub-levels with a damp- 
ing rate given by the relevant Larmor precession rate. 
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FIG. 8: Simulations of laser cooling CaF prior to Stark de- 
celeration in a travelling-trap decelerator. The points show 
the phase-space distributions in (a) the longitudinal direction 
and (b) the transverse direction. Green points: molecules en- 
tering the decelerator at ^ = 103 mm with no laser cooling 
applied. Red points: molecules entering the decelerator at 
z — 170 mm with laser cooling applied. Line: Phase-space 
acceptance of the decelerator for CaF(A'' = 1) decelerated at 
a rate of 2 x 10* m/s^. 



From these calculations we find that when the intensity 
at each of the 8 frequencies is /, and the frequencies have 
a common detuning, the photon scattering rate is like 
that of a two-level system with an effective decay rate 
roff/(27r) = 2.4 MHz and an effective saturation parame- 
ter Soff — ///s.cff, where /s,cff — 50mW/cm^. We use this 
effective scattering rate to calculate the scattering forces 
on the molecule. 

We take a supersonic beam of CaF with a mean speed 
of 340 m/s and a translational temperature of 3K. The 
spatial distribution of the source is Gaussian in all di- 
rections with a standard deviation of 5 mm in the longi- 
tudinal direction and 2.5 mm in the transverse direction. 
The beam passes through a 3 mm diameter skimmer at 
z = 68 mm and then into the decelerator. We consider us- 
ing the travelling-trap decelerator demonstrated in Berlin 
[55] to decelerate CaF molecules in the = 1 state, these 
being the ones that are subjected to laser cooling. In the 
travelling-trap decelerator the molecules are confined in 
a three-dimensional potential well throughout the decel- 
eration process and the instabilities that arise in the WF 



8 



decelerator [35] are absent. We take the geometry to 
be exactly that of |35j and the peak-to-peak amphtude 
of the apphed waveform to be 16 kV. In the case where 
there is no laser cooling applied we place the decelerator 
a,t z = 103 mm, close to the exit of the skimmer. In the 
laser cooling case the entrance of the decelerator is moved 
downstream to z = 253 mm. Cooling is applied in both 
transverse directions in the 150 mm gap between the exit 
of the skimmer and the entrance of the decelerator using 
retro-reflected laser beams with a common red-detuning 
of 7.5 MHz and with s^s = 0.5. Longitudinal cooling is 
applied using a laser beam that propagates along — z, has 
Seff = 0.5, and is pulsed on for 550 ^s starting 200 /iS af- 
ter the molecules are produced, so that they interact with 
the laser in the region between the skimmer entrance and 
the decelerator entrance. When this light turns on, the 
common detuning is -578 MHz, and then the frequency 
is chirped linearly at a rate of 0.048 MHz//is so that the 
total frequency change is 26 MHz. The decelerator traps 
have an acceleration of —2 x 10^ m/s^. Their initial speed 
is chosen as 340 m/s for the no-cooling case and 333 m/s 
for the cooling case. 

The results of these simulations are shown in Fig.jSj 
where we plot the phase-space distributions of molecules 
entering the decelerator in the two cases, and compare 
these to the phase-space acceptance of the decelerator 
in both the longitudinal and transverse directions. In 
the longitudinal direction the laser slows down molecules 
with speeds in the 335-345 m/s range, substantially com- 
pressing the speed distribution into a narrow bunch cen- 
tred near 335 m/s. The laser cooling does nothing to com- 
press the spatial distribution which continues to expand 
so that at the decelerator entrance the laser-cooled bunch 
is longer than the uncooled bunch. The over-all effect of 
the longitudinal laser cooling is to increase the number of 
molecules within the decelerator acceptance by a factor 
of 3.1. The picture is similar for transverse cooling: the 
cooling is very effective at reducing the transverse veloc- 
ity spread, but the spatial spread increases because of the 
extra distance the molecules have to travel. The number 
of molecules overlapping with the transverse acceptance 
increases by a factor of 1.4. With the laser cooling ap- 
plied in all three dimensions the number of decelerated 
molecules is predicted to increase by a factor of 6. The de- 
celerated molecules occupy a far smaller volume of phase 
space when the cooling is applied, because in the trans- 
verse direction the velocity spread is compressed to an 
area far smaller than the acceptance area. The simula- 
tions show that the phase-space density of the molecules 
exiting the decelerator increases by a factor of about 2000 
when the laser cooling is applied. So we see that although 
the laser cooling results in only a modest increase in the 
number of decelerated molecules, it greatly increases the 
phase-space density of slow molecules. A further increase 
in phase-space density could be obtained by using a sec- 
ond short pulse of longitudinal cooling as the molecules 
exit the decelerator. 



VII. CONCLUSIONS 

We have reported the Stark deceleration of CaF 
molecules for the first time and have compared the 
performance of two types of decelerator for slowing 
these molecules. The AG decelerator is used to slow 
down molecules in the strong-field seeking ground state, 
whereas the WF decelerator works only for molecules in 
weak-field seeking states and has been applied here to 
CaF in the (TV, M) = (4, 0) state. 

Using a 100-stage WF decelerator at 20 kV and at a 
phase angle of 75°, the molecules have been decelerated 
from 340 m/s to 224 m/s. A 200-stage decelerator of the 
same design would be able to bring these molecules to 
rest. For this case, the relatively high phase angle re- 
sults in a small acceptance. This could be increased by 
using a longer decelerator at a lower phase angle. A 
travelling-trap decelerator of the type recently demon- 
strated in Berlin |3S] is also suitable for decelerating 
CaF and should increase the number of slow molecules 
obtained. The WF or travelling trap decelerator is 
also suitable for slowing down molecules of considerably 
higher mass than CaF, although as the mass increases 
higher-lying rotationally-excited states are needed for ef- 
ficient deceleration, and the required decelerator length 
increases. The required length could be reduced by the 
use of cryogenically-cooled sources [37] to provide slower, 
more intense beams. Cryogenic sources to date produce 
either continuous beams or very long pulses that are not 
well matched to the small size of the potential wells in a 
typical decelerator. However, there does not seem to be 
any barrier to making far shorter pulses. 

Using a 21-stage AG decelerator at 20 kV, ground 
state molecules have been decelerated from 433 m/s to 
399 m/s. We observe a smaller decelerated signal than ex- 
pected from numerical simulations, and in the case where 
all 21 stages are used the decelerated bunch is moving a 
little faster than expected. The advantage of the AG de- 
celerator is its ability to focus molecules in strong-field 
seeking states. For large molecules, these are the only 
states available There are several difficulties in de- 
celerating large molecules to low velocity. Firstly, the 
molecules must stay inside the stability region for AG 
focussing and that requires either the lens length or the 
focussing power to change as the molecules slow down. 
One approach is to use long lenses at the start of the de- 
celerator and short ones at the end, with the long lenses 
segmented into a series of shorter ones so as to maintain 
the same deceleration per unit length. This arrangement 
has been discussed in detail in the context of decelerating 
ground state YbF [TH]. An alternative is to build each 
lens from four electrodes instead of two [24] in which case 
the focussing (and defocussing) direction of each lens is 
set by the polarities of the electrodes. This provides a 
more versatile structure because the exact arrangement 
of lenses is no longer fixed by the layout of the machine. 
Even with these more sophisticated lens arrangements 
it becomes very difficult to maintain stability of AG fo- 
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cussing at the lowest speeds [TH]- A further difficulty for 
strong-field seeking molecules is that AG deceleration is 
always accompanied by excess defocussing, as discussed 
above. This can be minimized by tapering the ends of 
the electrodes more gradually, for example by replacing 
the hemispherical ends of the electrodes in the present 
design with ellipsoidal ends. Finally, experimental work 
on AG deceleration and focussing consistently reports 
more loss than expected from detailed numerical sim- 
ulations, and these additional losses are thought to be 
due to mechanical misalignments. Exceptionally good 
alignment is needed to make the decelerator perform as 
expected. A cryogenically-cooled source operated in the 
effusive regime offers an alternative, and simpler, method 
to obtain slow-moving beams of large, heavy molecules. 

For some molecules, CaF being a good example, laser 
cooling promises to increase the phase-space density 
enormously. Laser cooling can both cool and decelerate 
the molecules. For the latter, a large number of photons 
have to be scattered, requiring a particularly 'closed' vi- 



brational structure. An alternative, requiring far fewer 
spontaneous emissions, is to combine laser cooling with 
Stark deceleration. We have simulated a possible exper- 
iment where transverse and longitudinal laser cooling is 
applied to CaF prior to Stark deceleration. The cooling 
is very effective in compressing the velocity distribution, 
but does not compress the spatial distribution. Instead, 
the size of the beam at the decelerator entrance increases 
because the decelerator needs to be moved downstream to 
allow space for laser cooling. In our simulations, the laser 
cooling increases the number of decelerated molecules by 
a factor of about 6, but increases the phase-space density 
of slow molecules by about 2000. 
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